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There are typographical errors in Eqs. (20) and (23), and some incorrect entries 
Table (1). I thank Ryan Lang for pointing these out. 

Equation (20) should read as follows: 



MJ \ r 
/ r \ 2 / 2M\ 

+3c2 (m) I 1 "-) 



M(M - r)(2M 2 + QMr - 3r 2 ) + 3 ^ 



r 2 (2M-r) 2 2 & \r - 2M 



Equation (23) should be replaced by the following: 
k 2 = (1 - 2Cf [2 + 2C (y - 1) - y] X 



!^2C (6 - 3y + 3C(5y - 8)) + 4C 3 [13 - lly + C(3y - 2) + 2C 2 (1 + y 
+ 3(1 - 2C) 2 [2-y + 2C(y - 1)] log (1 - 2C) 



The corrected values for the Love numbers in Table (1) are given in the table below. 
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Table 1. Relativistic Love numbers k 2 
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For a variety of fully relativistic polytropic neutron star models we calculate the star's 
tidal Love number k 2 . Most realistic equations of state for neutron stars can be approximated 
as a polytrope with an effective index n ~ 0.5 — 1.0. The equilibrium stellar model is obtained 
by numerical integration of the Tolman-Oppenheimer- Volkhov equations. We calculate the 
linear I = 2 static perturbations to the Schwarzschild spacetime following the method of 
Thorne and Campolattaro. Combining the perturbed Einstein equations into a single second 
order differential equation for the perturbation to the metric coefficient g tt , and matching the 
exterior solution to the asymptotic expansion of the metric in the star's local asymptotic rest 
frame gives the Love number. Our results agree well with the Newtonian results in the weak 
field limit. The fully relativistic values differ from the Newtonian values by up to ~ 24%. 
The Love number is potentially measurable in gravitational wave signals from inspiralling 
binary neutron stars. 



2. Introduction and Motivation 



A key challenge of current astrophysical research is to obtain information about the 
equation of state (EoS) of the ultra-dense nuclear matter making up neutron stars (NSs). 
The observational constraints on the internal structure of N Ss are weak: the observed range 
of NS masses is M ~ 1.1 — 2.2M ( iLattimer &: Prakashll2007l ). and there is no current method 
to directly measure the radius. Some estima tes using data from X- ray spectroscopy exist, 
but those are highly model-dependent (e. g. I Webb fc Barretl (120071 )). Different theoretical 
models for the NS internal structure predict, for a neutron star of mass M ~ 1.4M , a central 
density in the range of p r <y 2 — 8 x 10 14 gcm~ 3 and a radius in the range of R ~ 7 — 16km 
( ILattimer fc Prakashll2007l ). Potential observations of pulsars rotating at frequencies above 
14 00Hz could be used to constrain the EoS if the pulsar's mass could also be measured (e. 
g. Izdunik et all J^OoJ). 



Direct and model-independent constraints on the EoS of NSs could be obtained from 
gravitational wave observations. Coalescing binary neutron stars are one of the mos t impor- 
tant sources for ground-based gravitational wave detectors ( ICutler &: ThorneHl993l ). LIGO 
observations have e stablished upper limits on the coalescence rate per comoving volume 
( lAbbott et al.l 120071 ). and at design sensitivity LIGO II is expected to detect inspirals at a 
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rate of ~ 2/day (IKalogera et al.l 120041 ) . 



In the early, low frequency part of the inspiral (/ < 100Hz, where / is the gravi- 
tational wave frequency), the waveform's phase evolution is dominated by the point-mass 
dynamics and finite-size effects are only a small correction. Toward the end of the inspiral 
the internal degrees of freedom of the bodies start to appreciably influence the signal, and 
there have been many investigations of how well the EoS can be constrained using the last 
orbits and merg er, including constraints from the gravitational wave energy spec- 



scvcra 



trum (IFaber et al.ll2002l ) and from the NS tidal disruption signal for NS-black hole binaries 



( Vallisneri 20021 ) . Several numerical simulations of the hydrodynamics of NS-NS mergers 



have studied the dependence of th e gravitational wave spectrum on the radius and EoS (see, 



e.g 



Baumgarte & Shapiro! (120031 ) and references therein). However, trying to extract EoS 
information from this late time regime presents several difficulties: (i) the highly complex 
behavior requires solving the full nonlinear equations of general relativity together with rela- 
tivists hydrodynamics; (ii) the signal depends on unknown quantities such as the spins and 
angular momentum distribution inside the stars, and (iii) the signals from the hydrodynamic 
merger are outside of LIGO's most sensitive band. 

During the early regime of the inspiral the signal is very clean and the influence of tidal 
effects is only a small correction to the waveform's phase. However, signal detection is based 
on matched filtering, i. e. integrating the measured waveform against theoretical templates, 
where the requirement on the templates is that the phasing remain accurate to ~ 1 cycle over 
the inspiral. If the accumulated phase shift due to the tidal corrections becomes of order unity 
or larger, it could corrupt the detection of NS-NS signals or alternatively, detecting a phase 
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of the internal structure on the gravitational wave phase in this early regime of the inspiral 
is characterized by a single parameter, namely the ratio A of the induced quadrupole to 
the perturbing tidal field. This ratio A is r elated to the star's tidal Lo ve number k 2 by 
k 2 = 3GXR~ 5 /2, where R is the star's radius. iFlanagan fc Hindererl (120071 ) have shown that 



for an inspiral of two non-spinning 1.4M Q NSs at a distance of 50 Mpc, LIGO II detectors 
will be able to constrain A to A ^ 2.01 x 10 37 gcm 2 s 2 with 90% confidence. This number 
is an upper limit on A in the case that no tidal phase shift is observed. The corresponding 
constraint on radius would be R ^ 13.6 km (15.3 km) fo r a n = 0.5 (n = 1.0) fully relativistic 
polytrope, for 1.4M NSs (IFlanagan &: Hindererl 120071 ). 



Because neutron stars are compact objects with strong internal gravity, their Love num- 
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bers could be very different from tho se for Newtonian stars that have been computed previ- 
ously, e. g. by lBrooker fc Olid (119551 ). 



Knowledge of Love number values could also be useful for comparing different numerical 
simulations of NS binary inspiral by focusing on models with the same masses and values of 
A. 



In Flanagan &: Hindered (120071 ). the I = 2 tidal Love numbers for fully relativistic neu- 
tron star models with polytropic pressure-density relation P = Kp 1+1 ^ n , where K and n 
are constants, were computed for the first time. The present paper will give details of this 
computation. Using polytropes allows us to explore a wide range of stellar models, since 
most realistic models can be reasonably approximated as a polytrope with an effective index 
in the range n ~ 0.5 — 1.0 lLattimer fc Prakashl (120071 ). Our prescription for computing A is 
valid for an arbitrary pressure-density relation and not restricted to polytropes. In Sec. [31 
we start by defining A in the fully relativistic context in terms of coefficients in an asymp- 
totic expansion of the metric in the star's local asymptotic rest frame and discuss the extent 
to which it is uniquely defined. In Sec. HJ we discuss our method of calculating A, which 
is based on static li nearized perturbations of the e quilibrium configuration in the Regge- 
Wheeler gauge as in iThorne fc Campolattaro Jl967r i. Section [5] contains the results of the 
numerical computations together with a discussion. Unless otherwise specified, we use units 
in which c = G = 1. 



Definition of the Love number 



Consider a static, spherically symmetric star of mass M placed in a static external 
quadrupolar tidal field The star will develop in response a quadrupole moment Qij 0. In 
the star's local asymptotic rest frame (asympt otically mass c entered Cartesian coordinates) 
at large r the metric coefficient g tt is given by (jThornelll998l ): 



9tt) 



M 3Qi 



r 



2r 3 



n x n> - -5 tJ 



+ -S ijX l x j + O (r 3 ) 



(1) 



1 The induced quadrupolar deformation of the star can be described in terms of the star's I = 2 mode 
eigenfunctions of oscillation. 
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where ri 1 = x l /r; this expansion defines SijU and Qij- In the Newtonian limit, is related 
to the density perturbation Sp by 



Qi 



d x8p(x) I XiX 



-r 2 6, 
3 



(2) 



and Sij is given in terms of the external gravitational potential $ ex t as 



d 2 $ c 



'j 



dx l dxi 



(3) 



We are interested in applications to fully relativistic stars, which requires going beyond 
Newtonian physics. In the strong field case, Eqs. (j2J) and (131) are no longer valid but the 
expansion of the metric (CEJ) still holds in the asymptotically flat region and serves to define 
the moments Qy and £y. 

We briefly review here the extent to which these moments are uniquely defined since 
there are considerable coordinate ambiguities in performing asymptotic expansions of the 
metric. For an isolated body in a static situation these mome nts are uniquel y defi ned: £,;,- 
and Q jj are the coordinate independent moments defined by iGerochl (119701 ) and lHansen 
(11974J ) for stationary, asymptotically flat spacetimes in terms of certain combinations of the 
derivatives of the norm and twist of the timelike Killing vector at spatial infinity. In the case 
of an isolated object in a dynamical situation, there are ambiguities r elated to gr avitational 
radiation, for example angular momentum is not uniquely defined (iWaldl 1 1984 ). For the 
application to the adiabatic part of a NS binary inspiral, we are interested in the case of a 
non-isolated object in a quasi-static situation. In this case there are s till ambiguities (related 
to the choice of coordinates) but their magnitudes can be estimated (IThorne fc Hartldll985l ) 
and are at a high post-Newtonian order and therefore can be neglected. We are also interested 
in (i) working to linear order in 8^ and (ii) in the limit where the source of is very far 
away. In this limit the ambiguities disappear. 



To linear order in S^, the induced quadrupole will be of the form 



Qi 



-XSi 



(4) 



Her e A is a constant which is r elated to the I = 2 tidal Love number (apsidal constant) ki 
by ([Flanagan fc Hindered 120071 ) 

k 2 = ^GXR- 5 . (5) 



2 The I = 2 tidal moment can be related to a component of the Rieman n tensor R n n~,& of th e external 



pieces of the metric in Fermi normal coordinates at r = as — i?oiOj (see iMisner et al.l (|1973l )) 
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Note the difference in terminology: in Flanagan fe Hindererl (120071 ). A was called the Love 
number, whereas in this paper, we reserve that name for the dimensionless quantity hi- 



The tensor multipole moments and £y can be decomposed as 

2 



2m 
ij > 



m=-2 



and 



Qij ^ ^ Qm3^ 



2m 
ij > 



m=-2 



where the symmetric traceless tensors yf™ are defined by (jThornelll980l ) 
with n = (sin 9 cosy?, sin 9 sin (p, cos 9). Thus, the relation (j4]) can be written as 

Qm A£" m . 



(6) 



(7) 



(8) 



(9) 



Without loss of generality, we can assume that only one £ m is nonvanishing, this is sufficient 
to compute A. 



Calculation of the Love number 



4.1. Equilibrium configuration 

The geometry o f spacetime of a spherical, static star can be described by the line element 
jMisner et al.lll973h 

ds 2 = g ( $dx a dxP = -e u{r) dt 2 + e x{r) dr 2 + r 2 (d9 2 + sin 2 9d<j) 2 ) . (10) 

The star's stress-energy tensor is given by 

T a p = (p + p) u a u p + pgfj , (11) 

where u = e~ u ^ 2 d t is the fluid's four-velocity and p and p are the density and pressure. 



Nume rical integration of the Tolman-Oppenheimer- Volkhov equations (see e.g. iMisner et al. 
( 119731 )) for neutron star models with a polytropic pressure-density relation 



Kp 



l+l/n 



(12) 



where K is a constant and n is the polytropic index, gives the equilibrium stellar model with 
radius R and total mass M = m(R). 
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4.2. Static linearized perturbations due to an external tidal field 



We examine the behavior of the equilibrium configuration unde r linearized perturbations 
due to an external quadrupolar tidal field following the method of iThorne fc Campolattaro 
( I1967I ). The full metric of the spacetime is given by 



„ _ Jo) , . 

9 a/3 — 9 a/3 + 



(13) 



where h a p is a linearized metric perturbation. We analyze the angul a r dep endence of the 
components of h a p into spherical harmonics as in iRegge fc Wheeler! (119571 ). We restrict 
our analysis to the I = 2, static, even-parity perturbations in the Reg ge- Wheeler gauge 
(|Regge fc Wheeled! 19571). With these specializations, h a p can be written as (IRegge fc Wheeler 



1957 



Thorne fc Campolattardll967l ): 



^ = diag[e- , Wfl (r), e A(r) # 2 (r), r 2 K(r), r 2 sin 2 9K(r)] Y 2m (6, cp). 



(14) 



The nonvanishing components of the perturbations of the stress-energy tensor (fTTT) are 5T® = 
—5p = —(dp/dp)~ 1 5p and ST- = 5p. We insert this and the metric metric perturbation (1141) 
into the the linearized Einstein equation 5G^ = 87tST^ and combine various components. 
From 5G e e - 5G^ = it follows that that H 2 = H = H, then 5G r e = relates K' to H, and 
after using 5G e e + 5G^ = 16ir5p to eliminate 8p, we finally subtract the r — r component of 
the Einstein equation from the t — t component to obtain the following differential equation 
for H = H (for I = 2): 



H" + H' 



, ,'2m(r) t , . 
' ^r 1 + 47rr (p - p) 



+H 



6e> 



+ 4vre A 5p + 9p + 



P + P \ 
(dp/ dp) J 



a 



0. 



(15) 



where the prime denotes d/dr. The boundary conditions for Eq. (1151) can be obtained as 
follows. Requiring regularity of H at r = and solving for H near r = yields 



H(r) 



a r 



2tt 



l--^(0)+9p(0) +((W(W(0); 



?(0)+P(0)\ r2 + o , 3) 



(16) 



where a is a constant. To single out a unique solution from this one-parameter family of 
solutions parameterized by a , we use the continuity of H(r) and its derivative across r = R. 
Outside the star, Eq. ([TBI) reduces to 



h"+(--\')h'-(¥. + \*\h o. 

r J \ r 



(17) 
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and changing variables to x = (r/M — 1) as in iThorne fc Campolattarol (119671 ) transforms 
Eq. (1170 to a form of the associated Legendre equation with I = m = 2: 



1) H" + 2xH>-(q + -^—^H = Q. 



The general solution to Eq. (TIB"]) in terms of the associated Legendre functions Q2 2 (^) arid 
P 2 2 {x) is given by 



H = Cl Q 2 



21 r - l)+c 2 P 2 2 (^-l 



(19) 



M / \M 

where c\ and C2 are coefficients to be determin ed. Substituting the expressions for Q 2 2 (x 
and P 2 2 {x) from lAbramowitz fc Stegunl (11964] ) yields for the exterior solution 



H 



r \ 2 / 2M 



, r \2 / 2M 



M(M - r) (2M 2 + 6Mr - 3r 2 ) 3 

r 2 (2M-r) 2 + 2 ° k \ r - 2.U 



(20) 



The asymptotic behavior of the solution (1201) at large r is 



8 / M 



5 ^ r J 1 I \ r / 



3 <^ + °((f 



(21) 



where the coefficients c\ and C2 are determined by matching the asymptotic solution ( I2T1) to 
the expansion ([T]) and using Eq. ([H]): 



(22) 



We now solve for A in terms of H and its derivative at the star's surface r = R using Eqs. 
( 122|) and (120]) . and use the relation ([5]) to obtain the expression: 



A', 



8C £ 



1-2C) 2 [2 + 2C (<,-!)-</] x 



(23) 



!^2C (6 - 3y + 3C(5y - 8)) + AC 3 [13 - lly + C(3y - 2) + 2C 2 (1 + y 



+ 3(1 - 2C) 2 [2-y + 2C{y - 1)] log (1 - 2C) 



where we have defined the star's compactness parameter C = M/R and the quantity y = 
RH'(R)/H(R), which is obtained by integrating Eq. (1151) outwards in the region < r < R. 



-10- 



4.3. Newtonian limit 

The first term in the expansion of the expression (123]) in Mj R reproduces the Newtonian 
result: 

jv 1/2 -v 



k "-2{VTJ)' (24) 

where the superscript N denotes "Newtonian". In the Newtonian limit, the differential 
equation (ITS!) inside the star becomes 

r [dp/ dp r 2 ) ^ ^ 

For a polytropic index of n = 1, Eq. ( 125]) can be transformed to a Bessel equation with the 
solution that is regular at r = given by H — A^Jr/R J 5 / 2 (irr /R), where A is a constant. 
At r = R, we thus have y = RH'/H = (tt 2 - 9)/3, and from Eq. ([23]) it follows that 

k 2(n = l)=(-\ + ^) « 0.25991, (26) 



which agrees with the known result of iBrooker fc Olid (jl955l ). 



5. Results and Discussion 

The range of dimensionless Love numbers k 2 obtained by numerical integration of Eq. 
is shown in Fig. [I] as a function of Mj R and n for a variety of different neutron star 
models, and representative values are given in Table [2] These values can be approximated 
to an accuracy of ~ 6% in the range 0.5 < n < 1.0 and 0.1 < (M/R) < 0.24 by the fitting 
formula 

3 / 0.56 \ /M x ~ nm;; 



fe-2( _0 ' 41 + ^JUJ ' (27) 

Both Fig. [T]and Table [2] illustrate that the dimensionless Love numbers k 2 depend more 
strongly on the polytropic index n than on the compactness C = M/R. § This is expected 
since the weak field, Newtonian values k 2 given by Eq. (124]) just depend on n (through the 
dependence on y) . The additional dependence on the compactness for the Love numbers k 2 
in Eq. ( [23]) is a relativistic co rrection to this. For M l R ~ 10 -5 our results for k 2 agree well 



with the Newtonian results of IBrooker fc Olid (119551 ) . Figure [2] shows the percent difference 



3 Note, however, that LIGO measurements will yield the combination k 2 R 5 and therefore will be more 
sensitive to the compactness than the polytropic index. 
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(&2 — ^2)1^2 between the relativistic and Newtonian dimensionless Love numbers. As can 
be seen from the figure, the relativistic values are lower than the Newtonian ones for higher 
values of n. This can be explained by the fact that the Love number encodes information 
about the degree of central condensation of the star. Stars with a higher the polytropic 
index n are more centrally condensed and therefore have a smaller response to a tidal field, 
resulting in a smaller Love number. 

Some estimates of the masses and radii of neutron stars, given in Table [31 have been 
inferred from X-ray observations (jOzell 120061 ; I Webb fc Barretl 120071 ) using the information 
from three measured quantities: the Eddington luminosity, the surface redshift of spectral 
lines, and the quiescent X-ray flux. The range of the numbers A for these stars is shown 
in Fig. [31 LIGO II detectors will be able to establish a 90% confidence upper limit of 
A ^ 2.01 x 10 37 gcm 2 s 2 for an inspiral of two nonspi nning 1.4M W NSs at a dist ance of 50 



Mpc in the case that no tidal phase shift is observed ([Flanagan fc Hindered 120071 ) . 



The author thanks Eanna Flanagan for valuable discussions and comments. 
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Fig. 3. — The range of Love numbers for the estimated NS parameters from X-ray obser- 
vations. Top to bottom sheets: EXO0748-676, uCen, M 13, NGC 2808. For an inspiral of 
two 1.4M Q NSs at a distance of 50 Mpc, LI GO II detectors will be abl e to constrain A to 
A ^ 20.1 x 10 36 gcmV with 90% confidence iFlanagan k Hindererl 120071 ) . 
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Table 2. Relativistic Love numbers k 2 



n 


MIR 
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^2 
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0.5 


0.15 


0.173 
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Table 3: Estimated neutron star parameters from X-ray observations 



Cluster / object 


M(M ) 


R(km) 


M/R 


uj Cen a 


1.61 ±0.15 


10.99 ±0.71 


0.18 ±0.04 


M 13 a 


1.36 ±0.04 


9.89 ±0.08 


0.2 


NGC 2808 a 


0.84 ±0.12 


7.34 ±0.96 


0.22 ±0.01 


EXO 0748-676 b 


> 2.1 ±0.28 


> 13.8 ± 1.8 


0.2256 



"The parameters for these stars are the averages from the best fit values of the data in I Webb &: Barred (|2007l ) 
for their three different spectral fits. The errors given here reflect only the deviations among the best fit 
values for the fits. 

6 The values are taken ffom lOzell (|2006l ). 



Note. — Parameters used to generate Fig. ((3]). 



